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Abstract 

We prove smoothness of isometric immersions of surfaces en¬ 
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1 Introduction 

For a local isometric immersion u of a two-dimensional Riemannian man¬ 
ifold with positive Gauss curvature into there is a link between the reg¬ 
ularity of the metric and the regularity of u; in particular, if the metric is 
smooth then so is u. Without a priori assumptions on the regularity of u 
this link is broken; as shown by classical results of Nash and Kuiper even 
regularity is not enough. 
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In the present paper, we show that square integrability of the second fun¬ 
damental form of u is sufficient for the link to persist. In particular, if the 
metric is smooth, then u is smooth in the interior, provided that initially it 
belongs to the Sobolev space H^. 

We note in passing that replacing the assumption of positive Gauss curva¬ 
ture by zero Gauss curvature also breaks the regularity link. Indeed, local 
isometric immersions of class of the standard flat metric in into 
may fail to be see |Horllb] and |Horlla] for examples in this direction. 
Our regularity results for metrics with positive Gauss curvature rely upon 
earlier work by Sverak on the Monge-Ampere equation. Due to the low regu¬ 
larity, the passage from the scalar problem to the vectorial problem addressed 
here is not trivial. 

Relaxing regularity to regularity on the Sobolev scale is important for 
variational problems: the isometric immersions studied here arise nat¬ 
urally in thin him elasticity. In the present paper, we use this regularity 
result to derive homogenized bending models for convex shells from three 
dimensional nonlinear elasticity. 

Much work has been done on shell theories in elasticity. An overview of 
the derivation of models for linear and nonlinear shells by the method of 
formal asymptotic expansions can be found in |Giann] . In the case of linearly 
elastic shells, the models thus obtained can also be justihed by a rigorous 
convergence result, starting from three dimensional linearized elasticity. 

In the last two decades, rigorous justihcations of nonlinear models for rods, 
curved rods, plates and shells were obtained by means of F-convergence, 
starting from three dimensional nonlinear elasticity. The hrst papers in that 
direction are |ABP911 ILDR951 ILDR96j for strings and membranes (plates 
and shells, respectively). For plates, the nonlinear bending theory the Foppl- 
von Karman theory were derived in |FJM02j and |FJM06j . For shells, the 
nonlinear bending theory was derived in |F,lMMn^ and the von Karman 
theory in |LMPinj : see also |LMP11] for limiting models in an intermediate 
energy scaling regime between bending and von Karman theories, for the 
particular case of convex surfaces. 

Here we are interested in the derivation of the homogenized bending theory 
of shells by means of simultaneous homogenization and dimension reduction. 
Our starting point is the energy functional of three dimensional nonlinear 
elasticity: the elastic energy stored in a deformation u G of a 

reference shell C of thickness h > 0 around a surface S' C is given 


by 



( 1 ) 
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Here, is a non-degenerate stored energy function that oscillates periodi¬ 
cally in X, with some period £ 1. We are interested in the effective behavior 

when both the thickness h and the period £ are small. 

The separate limits h —)• 0 and £ —)■ 0 are reasonably well understood: in 
jFJMMO^ it is shown that, when H4 does not depend on £, then the func¬ 
tionals ([I]) T-converge, as h —)■ 0, to a two-dimensional bending shell theory. 
Regarding the limit e: —)■ 0, which is related to homogenization, the first rigor¬ 
ous results relevant in nonlinear elasticity were obtained by Braides |Bra85] 
and independently by Miiller |Miil87] . They proved that, under suitable 
growth assumptions on We, the energy ([T]) T-converges as e —)■ 0 (for fixed 
h) to the functional obtained by replacing We in ([1]) with the homogenized 
energy density given by an inhnite-cell homogenization formula. 

In this paper we study the asymptotic behavior when both the thickness h 
and the period e tend to zero simultaneously. Such a combination of di¬ 
mensional reduction and homogenization has been the subject of numerous 
papers: in the paper |BFF00j the authors study the same effects for nonlinear 
systems (membrane plate) by means of F-convergence, also without periodic¬ 
ity assumptions. In their general approach they also consider the possibility 
of oscillating boundary. The Id case of rods was addressed in |Nenl2j : more 
recently the plate model in the von Karman regime (see |NV13j l and in the 
bending regime (see |HNV14( IVellhj ) were analyzed. In these cases one does 
not obtain an inhnite-cell homogenization formula as in the membrane case 
(see [BFFOOj ). This is because the energy is essentially convex for small 
strain. Therefore one can use two-scale convergence techniques in all these 
cases. 

However, each case has its own peculiarities. In the von Karman theory of 
plates, one obtains a limiting quadratic energy density which is continuous 
in the asymptotic ratio 7 between h and e, for all 7 G [0, cxo]. Moreover, 
the case 7 = 0 corresponds to the situation when the dimensional reduction 
dominates and the resulting model is just the homogenized plate model in 
von Karman regime. The situation j = 00 corresponds to the case when 
homogenization dominates and the resulting model is the plate model in 
the von Karman regime corresponding to the homogenized functional. The 
same phenomenology is observed in bending regime for rods in |Neul2j . The 
limiting model of the homogenized plate in von Karman regime without pe¬ 
riodicity assumption is obtained in |Vel] and the model of rod in bending 
regime in |MVj . The nonlinear bending theory of plates is more involved in 
the periodic case. In |HNV14j the authors obtained asymptotic models in 
the case 7 G (0, cxo]. In | Veil 5] the author obtained the asymptotic model 
corresponding to the regime 7 = 0 under the additional assumption that 
-C h < 
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In the critical case h one may expect that the oscillations needed for 
the relaxation would have to include some other frequencies in addition to 
the one on which the material is oscillating. So two-scale analysis might not 
be sufficient in this case. In |N015] the authors derived the homogenized 
bending plate model by performing homogenization of 2D bending plate. 
The homogenization of shells in the von Karman regime (cf. [HV14j l also 
displayed a separation of the models for the case 7 = 0 . For generic shells, 
we were able to identify the model for the situations <C h -C £ and, 
surprisingly, also for the siuation h ~ Moreover, for convex shells we 
were able to completely resolve the case 7 = 0 . The geometric framework 
developed in |HV14j will be used in the present paper as well. 

As it was already emphasized in |HV 14j the bending theory for shells seems 
likely to be more involved. In this paper we address that case, but only for 
convex shells. To answer the question of homogenization of bending shell 
we need to answer on the question of additional regularity of isometries on 
convex surfaces, for which we apriori know (by the compactness result) that 
are FTThe main results of this paper are Theorem 12.II and Theorem l3.4[ We 
emphasize the fact that this regularity result is also needed for the simpler 
case of studying the shells that are layered in the thickness direction (this 
model also incorporated in the analysis here). For the plate model of layered 
materials in the bending regime, see [SchOTj). 

Let us explain the technical relevance of the regularity result Theorem 12.11 
for the homogenization. In the compactness result (Proposition 13.7p we are 
unable to identify — i.e., to obtain any additional information on — the 
matrix field B which arises as (part of) the weak limit of approximate 
strains. In the construction of the recovery sequence it is natural to try to 
recover this field by adding to the standard recovery sequence an additional 
term of order h; this term is called w in our proof of the upper bound. Under 
the convexity assumption, we can prove that with this term we can indeed 
recover any matrix held B. This completes the proof. 

In order to carry out this plan, one must analyze a system of equations 
whose coefficients depend on the limiting deformation, which is an isom¬ 
etry. The analysis of this system requires regularity of its coefficients, see 
Lemma 13.101 

No such regularity result was needed in (FJMMOS]. Instead, the gradient 
truncation method used in |F,IMn2j was employed. However, such a trunca¬ 
tion (Lusin type approximation) cannot be used in the present paper. 

Note that in contrast to the situation encountered here, in the case of plates 
(see |Sch071lHNV14] ) not every matrix held B can be recovered by adding 
an additional term. But one can recover any matrix held that is zero on all 
planar parts of the deformation. And this class of helds contains the ones 
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relevant for the relaxation. 


Notation 

Ck,a (denotes the space of fnnctions that are fc-times differentiable and their 
fc-th derivatives are Holder continnons with exponent a. W^’^, denote 
the standard Sobolev spaces. For H C M" we denote by Cq{A) the space 
of fc-differnetiable fnnctions with continnons derivatives that have compact 
snpport in A. We set Y = [0,1)^ and we denote by y the Enclidean space 
eqnipped with the torns topology, that is for all z E 1? the points y+z 
and y are identihed in 3^. We write C'°(y) to denote the space of continnons 
fnnctions / : R^ —)■ R satisfying f{y + z) = f{y) for all z ElA. We denote 
by C^iy) those fnnctions in C'*^(R^) fl C^iy) whose derivatives np to the 
k-t\i order belong to C°(y). We denote by H^{y) and H^{Sxy) 

the Banach spaces obtained as the closnre of C°°{y) and 
with respect to the norm in H^{Y) and H^{SxY)^ respectively. By 

L?‘{y),H^{y) etc. we denote the snbspaces of L?‘{y),H^y) etc. whose 
mid-valne over y is zero. For H C R^ measnrable and X a Banach space, 
L‘^{A,X) is nnderstood in the sense of Bochner. We identify the spaces 
L‘^{A, L‘^{B)) and L‘^{A x B) in usnal way. Standard basis vectors in R^ are 
denoted by e*. By SO(3) we denote the set of rotational matrices in R^^^, 
by so(3) the space of skew symmetric matrices in R^^^, while by Rgy^ the 
space of symmetric matrices in R^^^. For a matrix M by symM we denote 
the symmetic part of M, symM = \{M + M*). For A,B(Z R”, with A (& B 
we mean that A is compactly contained in B. 


2 Regularity of intrinsically convex sur¬ 
faces 


The pnrpose of this chapter is to prove the following resnlt: 

Theorem 2.1. Let 17 C R^ 6e a bounded domain and let g G C°°(17, R^yj^) 
be a smooth Riemannian metric on U. Assume that the Gauss curvature K 
of g is positive and let u G Hg{U), i.e., u G 17 ^( 17 , R^) solves 

(Vm)'^(Vm) = g almost everywhere on U. 


Then u G C°°{U). 

We use idea s an d a key result from the unpublished (but widely circulated) 


manuscript |Sve] . For our purposes, its main result is to deduce convexity of 
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solutions / of the Monge-Ampere inequality det V^/ > c > 0, cf. Lemma 
EH below. 

In Sve , this result is combined with a local graphical representation to 


prove smoothness of isometric immersions of subdomains of the sphere, 
endowed with the standard metric. 

Our proof of Theorem 12.11 also uses this idea of representing u locally as a 
graph of a function /. However, a priori u is not C^. Instead, we show that 
the normal to u is continuous. This is enough for u to be locally a graph. 
Finally, a bootstrap argument, using classical facts about two-dimensional 
Monge-Ampere equations on one hand and exploiting the link between u and 
its graphical representation on the other hand, implies that u is smooth. 


2.1 Continuity of the normal 

The purpose of this section is to provide a fairly self-contained proof of Propo¬ 
sition 12.31 below. In doing so, we combine ideas from [BN951 IBN96j and oth¬ 
ers, and we introduce a suitable notion of topological degree. For the reader’s 
convenience, we include proofs of its relevant properties. 

In what follows, we use the notation \g\ = det^f and = (detf/)^/^. The 
Christoffel symbols of g are denoted by F^^. Dehne h : B ^ ®sym by 
h = riu ■ V^u. The Gauss equation is easily seen to remain true for u G H^. 
It reads: 

dadfiU -f" T' (2) 

Since we are dealing with maps, we should verify the validity of Gauss’ 
Theorema Egregium. 

Lemma 2.2. If u & Hg{U) and h = ■ V'^u, then det h = K\g\ almost 

everywhere on U. 

Proof. As shown in [FJM061 Proof of Proposition 6], by approximation it is 
easy to see that the map u satishes 

\did2u\^ - didiu ■ d2d2U = ^ {d2d2gn + didig22 - 2did2gi2) (3) 

almost everywhere on U. Denote by P{x) the orthogonal projection from 
onto the subspace spanned by diu and d 2 U. Recall that h = Uu- V^u. Then 
we deduce from ([3]) that 

det h = —\P{did 2 u)\'^+P{didiu)-P{d 2 d 2 u)+- (cl 2 ^ 25 'ii + ^i^i 5'22 ~ 2 did 2 gi 2 ) ■ 

(4) 


6 











But in view of ([2]) we have 


P{d^df,u) ■ P{d,dsu) = 

We conclude that the right-hand side of (jl]) can be computed from g and 
its derivatives. Since g is smooth, the right-hand side of (jl]) agrees with 
K\g\. □ 

Proposition 2.3. Let U (ZM? he a hounded domain and let g G C°°{U) he 
a smooth Riemannian metric on U. Assume that the Gauss curvature K of 
g is positive and let u G Hg{U). Then the normal 

diu X d2U 
\diu X d2u\ 

to u is continuous on U. 


To prove Proposition 12.31 let i? > 0 and let (p G H^{B 2 R,Sf). Then, for 
almost every r G (0,i?) we have (p\dBr ^ hence by Sobolev 

embedding 

G C^idBr). (5) 

By a classical result of Schoen and Uhlenbeck, there exist (pk G C°°{B 2 r, S^) 
converging strongly in to ip. After possibly passing to a subsequence, we 
may assume that (p*, —)■ (p in H^{dBr) for almost every r G (0, R), hence 


(fk ^ (p uniformly on dBr- 


( 6 ) 


In fact, setting fk{r) = \Vp>k — V(pp, by the coarea formula we have 

|V(pfc - V(pp = / |V(pfc - V(pp0. 

Jbr 

Hence there is a subsequence such that fki^) —t 0 for Lf almost every r G 

( 0 ,/?). 

For the rest of this section we £x r G (0, R) such that (jl]) and ([H]) are satished 
simultaneously. Note in particular that ip{dBr) is compact, due to (jlD. Define 
the degree Q : —)■ M of with respect to B^ by setting 


ll/fclUpsK) = / dr 


'dBr 


Q{y) = f (V 

J Br 

where rj is any smooth 2-form on with Jg 2 h = 1 which is supported in the 
connected component Ay of \ ip{dBr) that contains y. 

We claim that d is well-dehned, i.e., that it is independent of the choice of 
g. We use the following lemma. 
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Lemma 2.4. Let A C 6e connected and let rj be a smooth 2-form on 
whose support is contained in A and which is such that = 0. Then there 
exists a smooth 1-form w on with support in A and such that rj = dw. 


Proof. The proof can be found in |Sch69] . 


□ 


In view of the lemma it remains to show that if ta is a smooth 1-form sup¬ 
ported in Ay then ip*{dw) = 0. 

Since dw is a 2-form and since ^ ip strongly in we see that 


(pl{dw) (p*{dw) strongly in L^{Br). 


Hence 



Due to (j6]), the compact set ipk{dBr) does not intersect the support of w for 
k large enough, because the latter has positive distance from the compact 
set (p{dBr). Therefore, (plw has compact support in Br- Hence, by Stokes’ 
theorem, the left-hand side of (jS]) is zero. This concludes the proof showing 
that Q is well-dehned by ([7]). 

Recall that the essential range oi ip\Br is the smallest closed set F such that 
<p>{x) G F for almost every x & Bp., as shown in [BN95] it is well-dehned. 
More or less directly from the dehnition of Q, we see the following: 

Lemma 2.5. Let (p G and define Q as in o. Then the following 


are true: 


(i) Q is constant on every connected component ofEf \ ip{dBr); 

(a) If Qiy) 7 ^ 0 then Ay is contained in the essential range F o/(p|b^. 

(Hi) Q takes integer values. 

Proof. To prove ([n]), assume that Ay is not contained in F. Then there exists 
an open set A <Z Ay\F and a normalized smooth 2-form rj supported on A. 
So (p*ri = 0 almost everywhere on B^.. Hence we would have Q{y) = 0. 

To prove (Imj) just note that the last convergence in ([H]) is also true for any 
other 2-form; in particular for the form 77 in ([7j). But for smooth (p, the 
right-hand side of ([7]), which is known to attain only integer values. □ 

Lemma 2.6. If p> = Uu, then we have Q > 0 on \ nu{dBr). Moreover, if 
Qiy) = 0 then Ay does not intersect the essential range of n\sBr- 






Proof. Denote by r]g 2 the standard area form on Then = Kd\o\g, 

due to Lemma [221 Applyig ([7]) with rj = prjs'^, we see that 



Q{y) 


{ponu)K dvolg 


' Br 


for every p G Since iL > 0 on 5^, we conclude that if Qipj) = 0 then 

p o n = 0 almost everywhere on Br- Since p was arbitrary, this implies that 
n[x) G \ Ay for almost every x G 3^. Since \ Ay is closed, by minimality 
of the essential range we conclude that it must be contained in \ Ay. □ 

Lemma 2.7. //n„ = p and if R > 0 is small enough, then Q is zero at some 
point in \ nu{dBr). 


Proof. There exists a constant C depending only on g such that 



K dvolg < CR^. 


We choose R so small that the right-hand side is bounded by 1/4 times the 
area of Ef. 


Assume for contradiction that Q 7 ^ 0 on E‘^\nu{dBr). Then by Lemma I2T] we 
know that Q is positive and so by Lemma 12751 we have Q > 1 on E‘^\nu{dBr). 
Since G H^{dBr), it maps dB^ into a set of zero area, cf. |Res89j . So the 
area of \ nu{dBr) is that of Hence there exist hnitely many pairwise 
disjoint connected components Ai, ...,Am of \ nu{dBr) and fji G C“(Aj) 
satisfying fji G [0,1], such that 

M 

5^Wy{V. = l})> (9) 

i=l 

Let Hi G Aj and note that ^ pointwise on Hence, recalling that 
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Q{.yi) > 1 , 



M „ 


/ {:ipionu) K dMo\ 
J Br 


K dvolg > E/ (^i o Uu) 


i=l 
M „ 




= 

i=l 


<(V'i^S2) 



M 



M 


M 




In view of ([H]) this contradicts our choice of R. 


□ 


Since Q vanishes at some point, it is in fact zero on a whole connected 
component A of the open set §>^\nu{dBr). But this means that (after possibly 
redehnig on a set of measure zero) does not take values in A. 

We assume without loss of generality that 63 G A and we denote by : 

\ {^ 3 } the stereographic projection. Since A C \ nu{dBr) is 

open, there exists p > 0 such that fl B 2 p{e^) does not intersect nu{Br). 



Hence \k o n G {H^ n M^). Since d' is 


conformal, we deduce from A' > 0 that the Jacobian of 'h o is positive 
(or negative) and bounded away from zero. Hence o is continuous, cf. 
|Res89] . (In fact, our arguments essentially even show that J/on^ is monotone 
in the Lebesgue sense, which quite directly implies continuity.) Hence Uu is 
continuous as well. 


2.2 Proof of Theorem 12.1 


Fix a point xq G H and let us prove that u is smooth in a neighbourhood 
of this point. We assume without loss of generality that xq = 0 and that 



where V = 63 ■ u and J/ : H —)■ is the 


in-plane component. 

Lemma 2.8. There exists a positive constant c such that det Vft' > c near 


0 . 
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Proof. We compute 


det V\I/ = 63 ■ 



X 



= 63 • diu X d 2 U = ^/g{nu ■ 63). 


By continuity of Uu and since > c, there is a neighbourhood of 0 on which 
the right-hand side is bounded from below by a positive constant. □ 

Since u is Lipschitz by hypothesis, we see that is Lipschitz and, by Lemma 
I2.81 that its Jacobian is bounded from below by a positive constant. Hence 
Clarke’s inverse function theorem |Cla76] shows that J/ is locally bi-Lipschitz 
on some ball B centered at 0. 

Dehne <h = and dehne / = H o $, which is a map '^{B) —)■ M. Then 
u{B) = graph/|^(s). For x G 'L(iJ) we dehne 

G{z) = {z,f{z)), 


so that u = G(d^). We assume without loss of generality that J/(0) = 0. 


Clearly, / G W^’°°{B). Moreover, / G by the chain rule for the compo¬ 
sition of Sobolev functions with a bi-Lipschitz map <F. The usual chain rule 
applies. We claim that / satishes the prescribed Gauss curvature equation 

detVV = ^(<^>)-(l + |V/P)2 (10) 


almost everywhere on d'(iJ). In fact, denoting the metric induced by G by 
^=(VG)^(VG) = J + V/®V/, 


and the normal to G hj riu = have 


Uu = 


■V/,1)^ 


(1 + |V/|2)V2’ 


( 11 ) 


because det^f = 1 -|- |V/p. Moreover, dad^G = (0, 0, Therefore, 

h = nu ■ V^G satishes 

vv 


a + iv/p)'/2' 

Taking determinants in flT^ . we see that 


( 12 ) 


detvy = (l + |V/|^)detft. 


(13) 
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Using the chain rule, it is easy to verify that 

h = (V'h)^ h(\I/) (VH/) almost everywhere. (14) 

A similar relation applies to g and 'g. Therefore, using Lemma 12.21 we see 
that flT^ implies flTU]) . 


Observe that the right-hand side of flTOD is positive and bounded away from 
zero and inhnity. Hence the following lemma implies that / is a (locally) 
convex function. 


Lemma 2.9. Suppose f G H^{B) satisfies detV^/ > c > 0. Then f is 
either locally convex on B or locally concave on B. 


Proof. This is proven in |Sve| . As observed in Sve , the results in 


(w hich were conjectured in jSvej 1 indeed allow to relax the hypothesis 

in jSve] to the hypothesis used here. We refer to |Horonl ILMPlSj and 
the references therein for more on this. □ 


We may assume that / is convex, and then, e.g. by a result of |Hei 6 n] or 
by |Oaf9nbl[Cm] . we may assume that / is strictly convex on B. One can 
then deduce from flTU]) that f G C^, but we give a direct proof of this: 

Lemma 2.10. We have f G (^^(^(i?)). 

Proof. We know that H —)• is continuous. And so is $. Since 

riu = we see that is continuous. Upon scalar multiplication of flTT|l 

with 63 , we have 

(l + |V/p)-V 2 = n^($).e3. 

Since n„(<h) • 63 is continuous and strictly positive, we conclude that (1 -I- 
is continuous, hence so is V/, by ffTTjl . □ 

The following lemma is |NS851 Theorem 1’]; cf. also |Caf90a] . 

Lemma 2.11. Let 0 < m < M < 00 and let F : Br —)■ [m,M]. If ip E 
C%Br) is a convex Aleksandrov solution of 

det V'^cp = F, 

then there exists p > 1 such that (p G W^^f{BR). Moreover, p —)■ 00 as 
M/m —)■ 1. 
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Since F is continuous, its oscillation is as small as we please on provided 
R is small enough. Hence Lemma 12.111 implies that there exists i? > 0 and 
p> 2 such that / G 

Since / G also G G And flT^ implies that h G 

LP{Bji). But V'L G L°° because u G Therefore, flTT|) shows that 

h G Lp^Br)). 

Since T is smooth and Vm is bounded, we deduce from ([2]) that u G W^’P{^{Br)). 
In particular, by Morrey-Sobolev embedding, there exists a constant 5 > 0 
such that T G {^{Br)) and / G C^’^{Br). Since K is, in particular, Lip- 
schitz on B, we have F G C^’^. E.g. by the results in |Sch82] . we therefore 
deduce from det V‘^f = F that / G Hence (IT^ shows that h G C^’^. 

Hence h G by (fTTD . Thus ([2]) implies u G for some positive 5 in 
some neighbourhood of the origin. 

Since on some neighbourhood of the origin u G for every constant unit 
vector e the function n = e • « is a solution of the Darboux equation 

det (V^n + T • Vv) = K\g\ (l — g~^ : (Vn ® Vn)) , 

where T ■ Vn denotes a suitable contraction of indices. This equation is 
elliptic with respect to v. Since g~^, T and K are smooth, and since v G 
we conclude by standard theory that v G C°°, hence u is smooth in a 
neighbourhood of the origin. 

Remark 2.12. Notice that if we only assume that g G then 

we can conclude that u G M^), for some 5 > 0. If g G 

for k > 3 then we can conclude that u G by standard elliptic 

regularity (see, e.g., IGTOll Lemma 17.16]). 

3 Homogenization for shells 

3.1 Preliminaries 

Unless additional properties are stated, a; C is assumed to be a bounded 
domain with boundary of class C^. We set / = (—|, \) and Of' = u x {hi), 
and H = uxl. From now on, S denotes a compact connected oriented surface 
with boundary which is embedded in For convenience we assume that S 
is parametrized by a single chart: From now on, f G denotes an 

embedding with ^(cj) = S. The inverse of f is denoted by r : S' ^ a;, and we 
assume it to be (7^’" up to the boundary. 

As before hj g = (V^)^(V^) we denote the Riemannian metric on u induced 
by f. Its Christoffel symbols are denoted by T^^^. In what follows we recall 
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some standard notions and set the notation. All notions are discussed in 
detail in most basic textbooks on Riemannian geometry. By we 

denote the isometries of the surface S. Notice that 

u € HUS) e S 6 HHU), 

where u = u o We repeat the geometrical framework we used in [HV14] . 

• The volume element on S is denoted by dvol^. 

• The scalar product on a vector space V is denoted by (x, y)y, and we 
define x ■ y = {x, y)^k, for x,y eM}. 

• We denote by TR the tangent bundle over S, i.e., the collection of 
tangent spaces T^S with x E S. A basis of T^S is given by ri(x), T 2 {x), 
dehned by 

Ta{x) = X E S. 

We can regard T^S as a subspace of then we write a -t = (a, 

• By T*S we denote the cotangent bundle. A basis of T*S is given 
by (r^(x),r^(x)) dual to (ri,r 2 ). It is uniquely determined by the 
condition 

r'^iry) = 

where Say is the Kronecker symbol. Observe that r“ is more commonly 
denoted by dx", but we will not use that notation. We can identify 
T^S with T*S via r i—)■ (r, •). Via this identification we can identify 
r G T*S with the unique vector v E T^S C with the property that 
T = {v, •)]g 3 on TxS. Notice that 

Vr = ( ^^2 ) ■ (15) 


By n : S' —)■ we denote the unit normal to S, i.e., 

ri(x) X r2(x) „ ^ o 

nix) = -— -tor all x E b. 

|ri(x) X r2(x)| 

We define = n. The normal bundle of S is denoted by NS and 

by definition has fibers N^S given by the span of n(x). We denote by 

Ts{x) = I — n(x) ( 8 ) n(x) 

the orthogonal projection from onto T^S. We will frequently deal 
with vector helds V : S' —)■ on the surface. 
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• The tensor product bundles TS^TS etc. are defined fiberwise. If T^S 
is regarded as a subspace of then T*S T*S can be regarded as a 
subspace of 

The symmetric product 0 F of two vector spaces (or bundles) E and 
F by definition consists of elements of the form 

1 

aQb := - (a 06 + 60 a) 

with a E E and b E F. 

Sections of the bundle T*S © T*S are called quadratic forms on S. 
Sections B of T*S 0 T*S can be regarded as maps from S into 
via the embedding l defined by l{B) = B(Ts,Ts). (On the right-hand 
side and elsewhere we identify (M^)* 0 (R^)* with By definition, 

B{Ts,Ts) : S -E R^^^ takes the vector fields v,w : S ^ into the 
function x e-)■ B{x)(Ts{x)v{x),Ts{x)w{x)). 

• For any vector bundle E over S we denote by L‘^{S]E) the space of 
all L^-sections of E. The spaces H^{S]E) etc. are defined similarly. 
Explicitly, e.g. for E = TS we have 

L\S,TS) = {f E L\S,R^) : f{x) E T^S for a.e. a; G S}. 

• For any vector bundle E over S with fibers Ex, we denote by L‘^{y, E) 
the vector bundle over S with fibers L‘^{y, Ex). The bundles H^{y, E) 
etc. are defined similarly. For example, L^-sections of the bundle 
H^{y,TS) are given by 

L\S,H\y,TS)) = 

{Z E L‘^{S,H\y,R^)) : Z{x) E H\y,TxS) for a.e. a; G S}. 

• For a scalar function / : S' —)■ R its gradient field along S will be 
denoted by df, which is also the notation for the corresponding 1-form. 
In other words, 

df{x){T) = Vr/(a;) for all r G TxS. 

Here and elsewhere V denotes the usual gradient on R^ (or on R^) of 
the extension of /, and Vr/ = t ■ V/ = . We extend these 

definitions componentwise to maps into R^. 

• For functions / G T^(S, H‘^{y)) we define Hessy / to be the section of 
the bundle E^ {y-,T*S Q T*S) over S given by 

{Ressy f){x,y) = (Vj/)o^(a;, |/)r“(a;) 0 r^(a;), 

where (V^/)^^ = dy^dy^f. 
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• For V G L‘^{S, we define the section Defyn of the bnndle 

L‘^iy,T*S QT*S) by 

(Def;yn)(a;,|/) = {sjmVyv{x,y))a^T°‘{x) t^{x). 

Here and elsewhere Vy is gradient in y with respect to the variable y. 

• The Weingarten map S of S' is given by S = dn, i.e., 


S{x)t = (Vrn){x) for all x G S, r G T^M. 


We extend S to a linear map on TS © S' = by setting S = S o T 5 , 
i.e., we define S{x)n{x) = 0. 

3.1.1 The nearest point retraction 

The nearest point retraction of a tnbnlar neighbourhood of S onto S will be 
denoted by tt. Hence 


7r(x + tn{x)) = X whenever |t| is small enough. 


After rescaling the ambient space, we may assume that the curvature of S is 
as small as we please. In particular, we may assume without loss of generality 
that TT is well-defined on a domain containing the closure of {x + fn(x); x G 
S, —1 < t < 1}, and that \Id + tS(x)| G (1/2, 3/2) for all t G (—|, \) and all 
X G S'. 

For a subset A C S' and h G (0,1] we define = {x + t?7,(x) : x G S, —h/2 < 
t < h/2}. In particular, the shell is given by 


S'^ = {x + tn{x) : X ^ S, t ^ {—h/2, h/2)}. 
We introduce the function f : S'^ —)■ R by 

t{x) = (x — vr(x)) • n(x) for all x G S'^. 


(16) 


We extend all maps / : S' —?• R^ trivially from S' to S'^, simply by defining 


/(x) = /(7r(x)) for all x G S'^. 


(17) 


In particular, we extend r, Ts and S trivially to S'\ i.e., we have S(x) = 
S( 7 r(x)) and Ts{x) = Ts{'k{x)) and r(x) = r( 7 r(x)) for all x G Sh The 
following identity is easy to verify and is given in |HV 14j 



(18) 


( 19 ) 
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3.1.2 Thin films 

To deal with thin hlms, we introduce the map S : a; x M —>■ by setting 

2 ( 2 ', ^ 3 ) = ^{z') + Z 3 n{^{z')) for all z' ^ u and z^ G M. 

As in |F,lMn2] we will use the diffeomorphism 0^ : —)■ fl given by 

Q^izi, Z 2 , Z 3 ) = {zi, Z 2 , zs/h), and for a map y : fl ^ we introduce the 
scaled gradient VhV = {diy, d 2 y, ^d^y). The counterpart of 0^ on the shell 
is the diffeomorphism 0 ^ : 5^ —)■ 5^ given by 

0^(a;) = 7r{x) + ^^^n{x). (20) 

It is easy to see that 

0 ^ o H = S o 0 ^ on fl/j. ( 21 ) 

For given u : 5^ —)■ we dehne its pulled back version u : —)■ by 

u = MoS. We also dehne its rescaled version y : —>■ by y{Q^) = u on 5^ 

and we dehne the pulled back version y of this map hy y = y oS. Then it is 
easy to see that 

(Why) oQ^ = Vu on Qh. ( 22 ) 

We dehne the rescaled gradient VhV of y by the condition 

(Vhy) oQ>^ = Vu on (23) 

This corresponds to the fact that the rescaled gradient on th canonical do¬ 
main is just the physical gradient (on the physical domain) translated to the 
canonical domain. Using fl2l|l and fl2^ it is easy to see that 

V,^=Vu(S)((VS)o(0'^)-i). (24) 

Since Vf = u, formula ffT5]l and formula ffT^ show (recall that n is extended 
trivially to S^)\ 

V0^ = Vvr -I- —Vu + —n®n = Ts{I + tS)~^ + —ST 5 (/ -|- tS)~^ + —n ® n. 
n h n n 

Since Ts clearly commutes with S, we see that Ts commutes with (/ -|-tS)“^ 
as well. Hence 

V0'* = (4 + ^S)(/ + fS)-^ on (25) 

where Ih = Ts + j^n<^ n. To express VhV in terms of Vy, insert the dehnition 
of y into fl25]) and use fl25D to hnd 

VhV = Vy {h + tS){I+ thS)-^ on S\ (26) 
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3.1.3 Two-scale convergence 

Recall that we extend the chart r trivially from S to S^. We make the 
following dehnitions: 

• A sequence (/^) C is said to converge weakly two-scale on to 

the function / G (y)) as h —)• 0, provided that the sequence 

(/^) is bounded in and 

lim f 

J 

for allpe C'0(5\C'0(y)). 

• We say that strongly two-scale converges to / if, in addition, 

lim ||/^||l 2 ( 51 ) = ||/||L2(5ixy). 


(x) p{x, r{x)/e) dx = 



51 Jy 


f{x,y) p{x,y)dydx. (27) 


We write denote weak two-scale convergence and /^ -^ / to denote 

strong two-scale convergence. If ^ / then ^ Jy f{-,y)dy weakly in L^. 
If is bounded in L‘^{S^) then it has a subsequence which converges weakly 
two-scale to some / G L‘^{y)). These and other facts can be deduced 

from the corresponding results on planar domains (cf. |A1192l IVisObj i by 
means of the following simple observations. 

Dehning /^ = /^ o H and f{z,y) = f{E{z),y), and taking 

p(A2/) = p{^{z),y) {detVE'^{z)VE{z)y^‘^ , 
a change of variables shows that fl27D is equivalent to 

f f^{z)p{z,z'/e) dz ^ [ [ f{z,y)p{z,y) dy dz, (28) 

Jn JnJy 

where z' is the projection of 2 ; onto M^. Hence ^ / on if and only if 
y/* A j on 12 in the usual sense. 

When /^ : S' ^ M, then f^-^fonS means, by definition, that the trivial 
extensions converge weakly two-scale on S^. In particular, /^ -^ / on S' if 
and only if /^ -^ / on cn. All these dehnitions are extended componentwise 
to vector-valued maps. For quadratic forms q, on S we say g if 

q^ij^a) qij^cr) for all r, a G C^(S', TS'). A similar dehnition applies 
to other bundles. In the end we give the compactness lemma which gives 
the information about two-scale limit of scaled gradients. For the proof see 
|HV14l Lemma 4.3]. 
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Lemma 3.1. Let w’^ G be such that 

limsup (+ \\Vhw^Wl^^s^) ) < oo. 

h-i-O 


Then there exists a map wq G M^) and a field e L‘^{S x I xy, 

of the form 


w,eL\S;H\y,R^)) 
W2 G L‘^{S xyx 


= < 


{VyWi^w^) for some 

z/7 = 0, 

VyWi, -d^Wi ') for some Wi G L^(S'; H^{I x y, R^)) 


{VyWi,W2' 


if-f E (0,oo), 

wi G L\S X I,H\y,R^)) 
W2 G L‘^{S X /,R3)) 


for some 
i /7 = cx), 

such that, up to a subsequence, —)■ Wq in and 


(29) 


3 

VhW^ dwo oTs + {Hfijijx'' ® weakly two-scale on S^. 

i,j=l 


Here, wq is the weak limit in H^{S) of JjW^{x-\-tn{x))dt and G x 

j;.]^ 3 x 3 ) defined by H^{x,y) = H^{Ti{x),t{x),y). 


3.1.4 Elasticity framework 

Throughout this paper we assume that the limit 

7 := lim -fy 

exists in [0, cx)]. We will frequently write e instead of e{h), but always with 
the understanding that e depends on h. In the case 7 = 0 we will additionally 
assume that e{h)^ h. 

From now on we £x a Borel measurable energy density 

W : xR^ X R3^3 ^ [0, CX)] 

with the following properties: 


19 



(Wl) (continuity) W{-,y,F) is continuous for almost every y e and F G 

^3x3_ 

(W2) (periodicity) W{x,-,F) is 3^-periodic for all x G 5^ and almost every 
F G 


(W3) (objectivity) For all (x, y) G 5^ x 3^ we have hF(x,j/,/) = 0 and 
W{x, y, RF) = fF(x, y, F) for all F G M3x3^ ^ ^ go(3). 

(W4) (uniform quadratic growth near identity) There exist constants 0 < 
a < (3 and p > 0 such that for all (x, p) G 5^ x ^ we have 

fF(x, y,F)>a dist2(F, S0(3)) for all F G 

W{x, y,F)<l3 dist2(F, S0(3)) for all F G with dist^(F, S0(3)) < p. 


(W5) (uniform expansion)For each (x, p) G 5^ x 3^ there exists a quadratic 
form Q{x, y, •) : —)■ M such that 


ess sup 

(3;,j/)eSi xV 


\W{x,y,I + G)-Q{x,y,G)\ 

|G|2 


— y 0 as G — y 0. 


(30) 


In the following lemma we collect the properties of Q{-, ■, ■). They are direct 
consequence of the properties of W. For a proof see |Neul21 Lemma 2.7]. 

Lemma 3.2. Let W be as above and let Q be the quadratic form associated 
to W through the expansion fl30|) . Then 

(Ql) Q{-,y, •) is continuous for almost every y G 

(Q2) Q{x, •, G) is Y-periodic and measurable for all x E fl and all G G 

(Q3) for all X E fl and almost every y E the map Q{x,y, •) is quadratic 
and satisfies 

a|symGp < Q{x,y,G) = Q{x,y,sjmG) < (3\sjmG\‘^ for all G E 

The elastic energy per unit thickness of a deformation G of 

the shell is given by 

j'*(m'^) = 1 f iy(0^(x),r(x)/5,VM^(x)) dx. 
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In order to express the elastic energy in terms of the new variables, we asso¬ 
ciate with y : ^ the energy 

-^^( 2 /) = [ W {x,r{x)/e,Vhy{x))det{I+ t{x)S{x))~^ dx 

Js^ 

W {x + tn{x),r{x)/e, Vhy{,x tn{x))) dt dYo\s{x). 

By a change of variables we have 

= i [ W {x,r{x)/e,Vhy^{x)) Idet V(0'*)-^(a:)l dx, 
h Jsi 

where y^{E^) = u^. Using fl2^ it is easy to see that 

= jh^yh^ (1 + o{h)) as h ^ 0, (31) 

where \0{h)\ < Ch. 

3.2 Asymptotic energy functionals and main result 

Next we will introduce the asymptotic energy functionals. In order to do so, 
we need the dehnition of the relaxation helds and the cell formulas. 

Recall that a Q b = b + b <S) a). We now dehne relaxation operators 

with range in the space of L^-sections of the vector bundle over S with hbers 
given for each a; G S' by 

(/ X y, {T*S O T*S) © {T*S 0 N*S) © {N*S © N*S )). (32) 

Of course each of these hbers is isomorphic to L^(/ x (V; We now make 

the following dehnitions: 

Set D{Uo) = x H^{y) x L‘^{I x y,R3) and dehne 

^oiC^ 9 ) = Def;y C + 25 ( 0 r" © n + g^n © n — f Hessy (p for all {(, ip, g) G h^(S, D{Uq)). 

Set D (Woo) = X L‘^{I,H^{y)) x and dehne 

Woo(C, P-, c) = Defy C+ 2 ( 9 y,,p+CQ,)r"©n+C 3 n©n for all (C, p, c) G (S', D (Woo)) • 

Set D{U^) = H\I X y;M2) x H\I x y) and dehne 

^■yiC, p) = Defy C + {dy^P + -^3(0)7-“ © n + {-d‘ip)n © n 

7 7 
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for all (C,p) e L^{S,D{U^)). 

By trivially embedding D{Uq) as constant maps into D(Uo)), we can 

regard Uq also as a map from D(Uo) itself into ([32]). For each x E S the 
hberwise action Wg of Uq is 

= (DeiyC){x)+2gaT°‘{x)Qn{x)+g3n{x)(^n{x) for all {C,g) G D{Uo). 
For each x E S we dehne L^\l x 3^) = uI^\d{Uq)), i.e., 

Lt\i xy) = [ul,^\c,g) : (C, <?) e D{Uo)] . 

This is a subspace of fl3^ . i.e., of L^(/x3^; Mgym). We denote by Lq{I xy) the 

vector bundle over S with hbers L^q\i xy)] in what follows we will frequently 
omit the index (x) for the hbers. The bundles Ty(/ x (V) , for 7 G (0, cxd] are 
dehned analogously. The elements of these spaces are the relaxation fields. 
For 7 G [0, cx)] and x G S' we dehne Qy{x) : T*S ® T*S ^ R, hj setting 

Qy(x,g)= inf / / Q(x + tn{x),y,p + tq + U{t,y)')dydt, (33) 
ueLfihixy) Jl Jy ^ ^ 

peTiStSTiS 

for each x G S' and q E T*S ®T*S. For x E S and q E T*S QT*S dehne the 
homogeneous relaxation (cf. [LMPinj h 

Q 2 {x,t,q)= min {Q{x + tn{x),M) : M{Ts,Ts) = q{Ts,Ts)}. (34) 
Then it is easy to see that 

Qo{x,q) = inf / Q 2 [x + tn{x),y,p + tq + Defy ( - tUessy (p) dt dy, 
Jixy ^ ^ 

(35) 


where both inhma are taken over all C, E II^(y,R^), all p E H‘^{y) and all 
p G T*S © T*S. In the case when Q does not depend on t , i.e., the material 
is homogeneous in the thickness direction we hnd: 

Qo{x,q) = ^ inf / Q 2 {x,y,q+ }lessy p) dy. 

12 ipeH^iy) Jy 

Remark 3.3. The formulae fl34D and fl3TD eorrespond to the fact that in 
the regime when dimensional reduction dominates, one can, at least on the 
level of energies, first do the dimensional reduction fl5T|) and then the homog¬ 
enization dSSI)- However, doing homogenization of the bending shell model 
is a different problem, since one has to deal with the restriction of being an 
isometry (see \N01h^ for the planar case). 


22 






As in |NV13j . one can prove pointwise continuity of the energies in the pa¬ 
rameter 7 , i.e., that for all a; G S, q E T*S QT*S we have 

lim Q.y{x,q) = Qoo{x,q) and \miQJx,q) = Qo{x,q). 

■y—>-oo 7—>-0 


For u G Ffjfo) by we denote the relative change of Weingarten maps 

S:: = iRoE~yViRoE-^)Ts-n (36) 


= [Ro] 
= {Ro: 




:-l\t 


V((floH-‘)n)rs-S 

)‘S„,s|(iioS-‘)Ts-S, 


(37) 

(38) 


where R = Vu G H^{oj] S0(3)) and we have denoted by V{R^ o S ^)Ts ■ n 
the operator on given by 


{y{R^ 




)Ts ■ n)v = 'SJ{R^ o S jTsv n, for all v G 


Notice that G L\S]T*S QT*S). 

For 7 G [0, 00 ] we dehne the functionals : H^{S; M^) —)■ M by setting 


I,{u) = ^ 


if M G H? 


ISO’ 

otherwise. 


/5Q7 (x,S(;(a;))dvol5(x), 

-Foo, 

We are now ready to state the main result of this part. 

Theorem 3.4. Let W be as above and assume that G H^{S^;M^) satisfy 

—2 jhf^.h\ 


lim sup h ^J'^{u ) < 00 . 

h^O 


(39) 


Define ^ by y^{Q^) = u’^. Put 7 = hm/i_^o ^ <ind in the case 

7 = 0 assume additionally e{hfi h. Then the following are true: 

(i) (lower bound). If y^ is such that y^ — y^ —)■ u, strongly in then 
we have 

liminf > iRu). 

u •. f\ ^ ' I \ / 


(a) (recovery sequence). Assume additionally that oj is star shaped with 
respect to some ball B <Z u and that S is convex, i.e., there exists 
C > 0 such that for every x E S we have 

— |r|^ < S(x)r • r < C\t\^, Vr G T^S. 

o 

Then we have that for any u E there exist E satis¬ 

fying (gg), and such that y^ —?• u, stronly in H^. Moreover, we have 

lim h-^.J^{u^) = Uu). 

u ^r\ ^ ' / v / 
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This is complemented by the compactness result, given in |FJMM0^ Theo¬ 
rem 1]. 

Theorem 3.5. Suppose a sequence G has finite bending en¬ 
ergy, that is is satisfied. Define ^ as above. Then there 

exists u G such that on a subsequence 

y^ — y^ dx ^ u, strongly in ,R^), 

Js^ 

Vhy^^R, strongly in L‘^{S^ 

as h ^ 0 after passing to subsequences and extending u and n trivially to S^. 
Here R G H^{S, SO(3)) is defined by RTs = VuTs and Rn = Rti x Rt 2. 

3.3 Two scale compactness and lower bound 

We impose that the sequence {vfi) C satisfies 

limsup < oo, (40) 

h^O 

and for a given sequence (pR) C we define the sequence {y^) C 

M^) by the expression y^{Q^) = u^. The following lemma gives us the 
compactness result we need. It is a direct consequence of [FJM02( Theorem 
3.1] and of the arguments in (FJMQSj. We refer also to (FJMMQS]. The 
result is obtained by dividing the set oj into the small cuboids of size 

d{h) in the in-plane direction and height h and then applying the theorem 
on geometric rigidity on the domains S(C'f) for every i and h > 0 and by 
smoothening. 

Lemma 3.6. There exist constants C,c> 0, depending only on u, such that 
the following is true for every h > 0 sufficiently small: if u E 
then there exists a map R : oj ^ SO {3) which is piecewise constant on each 
cube X + 5{h)Y with x G and there exists R G such that 

for each ^ G which satisfy |^|oo = niax{|^ • ei|, |^ • 62 !} < h and for each 
OJ (Z u which satisfy dist(a;, du) > ch we have: 

II Vu O (0'^)-! O S - + ||i? - .Rlli.(^) + R\\R- 

+A=|| < C Ijdist {Vu o (e*)-' o H, SO(3)) 11^,, 


24 










where {t^R){z') := R{z' + 0- Here 


5{h) := < 


eih), if lim —— > 0 
^ ^ h^oe{h) 

h 

*/lim ^TT = oo, 




h — ^0 s(ytl^ 

otherwise. 


For a sequence {u^) that satisfies iOl), we define the sequence {G^ 

L\S^- 


C 


of approximate strains with 

(i?'* o o (O'*)-! - / 


= 


h 


(41) 


where R^ : oj ^ SO(3) is a sequence of piecewise constant maps given by 
Lemma l3.6[ Notice that o (0^)“^ = VhV^, where the operator is 
given by ([23]). 

For u (£ u we denote by S' = S(a;) and by S^ = a; x /. We are interested 
in identifying two-scale limits of the sequence (G^). We will not give the 
complete proof, only reduce the case to the planar situation , when the 
rest of the proof goes by not so difficult adaptation to the proofs given in 

|HNV141IWT^ . 

Proposition 3.7. Let {u^) he a sequence of deformations that satisfies m. 
let G^ be defined with (^7|) and let u G be determined by Theorem AS.f^ 

Then for every u (£ u such that duj is G^’^ we have the following convergence 
on a subsequence: 

• case 7 G (0,oo).- there exist B G L‘^{S,T*S © T*S) and (C,p) G 
L‘^{S, D{U.y)) such that 


symG'^ ^ B + tS[^+l{^{C,p)-, 


• casey = oo: there exist B G L'^{S,T* SqT* S) and {(, p, c) G Lf{S]'D{Uoo)) 
such that 

sjmG^ ^ B + tSl+Uoc{C,p,c), 

• case 7 = 0, -C h <C there exist B G L^(S';T*S' 0 T*S) and 

{C,Tt9) ^ D{Uq)) such that 

symG'* ^ B + tSl+Uo{C,T,9)- 
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Proof. Take u ^ u such that u has boundary We decompose the defor¬ 
mation on the set = u x I, where a; a; is such that du is 

y^{x + tn) = y^{x) -|- th{R^ o 'Er^)n + hz^{x tn), (42) 

where is given by Lemma ITHI and y^{x) = Jjy’^{x + tn{x))dt. Notice that 

J z’^{x + tn)dt = 0. (43) 

We compute using ffT^ and fl26|l 

V;,/ = Vy’^Tsil+ thS)-^+ {R’^oE-^)n(^n (44) 

+thV{R’^ o E-^)Ts{I + thS)-^ ■ n + th{R’^ o E-^)STs{I + thS)-^ 

+hVhZ^. 


(45) 


(46) 


Notice that 

Vf = Jvy’^Ts{I + tS)dt. 

From Lemma 00] and 03.61) we conclude that 

\\Vhy’^-R’^oE-^\\^,^^,^<Ch. 

From the identity fl26|) and 0451) we also conclude that 
\\Vf-~{R^oE-^)Ts\\^,^s^^<Ch. 

From the last two expressions and 0441) we can conclude that 

limsup \\VhZ^\\^ 2 ( 01 ) < oo. 

h^O 

By using Poincare inequality and 045]) we also conclude that 

limsup < cxD. 

h^O ^ 

In order to identify the two scale limit of the approximate strain we decom¬ 
pose it in the following form 

V^'*Ts(J + thS)-^ - {R'^ c 


= {R^o 

+{R^ 


-^)Ts 


.-utiR 


^ o 


h 

{R^ 




h 


- n (X) n 

+t{R^ o E-yV{R'^ o E-^)Ts{I + thS)-^ ■ n 
+t{R^ o E-^Y{R^ o E-Y^Ts{I + thS)-^ 

+ {R^ oE-^'V Hz\ 
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where, similarly as before, by V{R^ o H ^)Ts{I + thS) ^ ■ n we have denoted 
the operator on given by 

{V{R’^oE-^)Ts{I+thS)-^-n)v = V{R'^oE-^)Ts{I+thS)-\ n, for all v e Ml 


Notice that 
V^"Ts(J + thS)-i 
h 


{R’^oE-^)Ts 


+0(h^), 


where ||0(h^) 11 ^ 2 ( 51 ) < Ch'^. Notice that Vy^Ta = daV^^ for a = 1,2. Since 
now we can, by pnlling back, easily translate the problem to the planar one 
by using the techniques developed in |HNV14] (see the proof of Proposition 
3.2) and Lemma [3.II we can conclude the following 


• case 7 G (0, oo) 

There exist G w,v e L^S, H\y,R^)), B G 

z G X 3^,M^)) such that 


VtTs-{R'^oE-^)Ts 

h 


*<3;j<2 

+-(2/-V')i?oS-'(a;)Ts 

7 

+ dy.Vi{x,y)R 

j<3;i<2 


{R^ O S-p - (R^ O S-p 2 
- h -" - 

y{R^ oE-^)Ts{I+ thS)-^-n - 




3 

i=l 

+^{y ■V')RoE-\x)n; 

V{RoE-\x))Ts-n 

+ ^ dy^Wi{x,y)R 
*<3;i<2 
3 

( Vyz{x, t, y), ^dszix, t,y)^^R 

*j=i 


Here R is limit of the sequence i?^, i.e., R^ and we denoted by {y ■ 
VR) = yidiR + 2 / 2^2 -R- Using the fact that R^daR for a = 1, 2 is skew 
symmetric matrix we conclude that there exist B G L‘^{S,T*S QT*S) 
and (C,p) £ L'^i.S.DiJA^)) such that 


symG"-i? + tS:: + W,(C,p), 
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where 


B 


o,/3=l,2 

( {MB):,, \ 

M{z + v + tw^ + tw) + 7t f [MB),2 1 > 


and where M G L°°{S-, = {RoE ^)V-^ ■ t,. 


• case j = oo. In the same way as before we can conclnde that there 
exis^^o e L^(S,R^), w,v e L^(S, ff\y,R^)), B e L^S,R^^^), z E 
L^(S X /, M^)), dEL'^iSxI, R3) such that 


Vy’^Ts-{R^oE-\x))Ts 2 
h 


{R^ o S-i) - (R^ o S-i) 2 

- h -" - 

V(^'^oS-^)Ts(/ + thS)-^-n - 




Bij{x)R^T^ 

*<3;j<2 

+ dy.Vi{x,y)R (S)T^] 

*<3;i<2 

3 

i=l 

V{RoE-\x))Ts-n 
+ Y dy.Wi{x,y)R 

*<3;i<2 

3 

Y {^yZ{x,t,y),d{x,t))^.R (^T^. 

*J=1 


From this we conclnde that there exist B E L‘^{S-,T*S © T*S) and 
{C,p,c) E L?‘{S]V{Uao)) snch that 


where 


symG"-i? + tS:;+Woo(C,P,c), 


B 



c 


r“©r^ 

/ l{wi + vi) \ 

M I l{w2 + V2) + Mz, 

( \{di + w\) \ ( i(M5)3i \ 

M ( 1(^2 + ^2) ) ( \{MB )32 1 ) 



and M is defined above. 

The case 7 = 0 , s{hY -C h e{h) reqnires different approach. We will 
explain how to adapt the approach from |Vell5] . Firstly by nsing we 
conclnde that 


IIV'/ - o0,i?"(r2 o0)||l2(<,;m3x2) < Ch. 

Denote by = {R^{ti o ^), i?^(r 2 o ^)). We look the following minimization 
problem 

min / iW'v^ - dx'. 

i;'*eiri(s.R3) Jg 
J- 

The associated Enler-Lagrange eqnation reads 

j-A'v^= - V-p^ inD 

1 duv’^ = p^ ■ v on 

snbject to J~v^dx = 0. Above, z/ denotes the normal on du. We obtain 
by standard regnlarity estimates that G nnder the assnmption 

that du is C^’^. In this case it holds 

II'^^^I|h 2(S;;R3) < C'dI V ■ P^||l2(^.r3) + ||p^||2,2(^.K3)) 

< Ch, 

||V'F"-V'^"|U2(^;M3) < Ch, 
||V'F"-(.R"(rioO,^'^(r2oO)|U2(^;K3x2) < Ch, 

for some C > 0. We now decompose the seqnence of the deformations {y^) 
in the following way 

y^{x + tn) = v^{x) + {y^{x) — v^{x)) + th{R^ o E~^)n + hz^{x + tn), (47) 
where o We also decompose the seqnence of strains (G^) in the 
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following way 


= T Y1 {{R’"°^~^)ra-Vv^Ts{I + thS)-W0)T‘^^r^ 

«,/3=l,2 

+ ^ ^ ((i?'* o S-i)n • Vv^Tsil + thSyGo,) n®T^ 

+1 (R’^ O E-^y(yfTs{I + thS)~^ - Vv^Tsil + thS)-^) 
h ^ 

iR’^ oE-y\R’^ oE-y - I 

+ - -- - - n®n 

h 

+t{R^ o E-y^ViR'^ o E-yTs{I + thS)-^ ■ n 
+t{R^ o E-y\R^ o E-ySTs{I + thS)-^ 

+ {R^ oE-y^WhZ^. 

By using the tools developed in [VellSj for the planar case, adapted to this 
framework, we can conclude that there exist G L‘^{S,T*S QT*S), G 
L2(5, w G L2(5, H\y, M2)), pe L\Sxy,R^),k e L^S; H^y, M^)), 

r G L^{S, H\y, M^)), if e H‘^{y), cf) e L\S X /, H\y, m3)), deL^SxIx 
3^, m 3) such that 


i ^ {{R^ o H 3)ro • dr^v^) ®T^ ^ 

o,/3=l,2 

^ ^ {{R^ O E-^n ■ dr^) n®T^ ^ 
* 0 = 1,2 

^{R^oE-y\VtTs-Vv^Ts) - 

{R’^oE-y\R^oE-y-I 2 
h "" ^ 

V{R^ oE-yTs-n ^ 

Vhz^ ^ 


B^{x) + Def^; w, 

Pa{x,y)n® for a = 1,2, 
B'pxyisTi 

*<3;i<2 

+ ^ dy.ki{x,y)R ®Ty 

j<3;i<2 

3 

-RHess^; </?, 

3 

^ {Vy(j){x,y),d{x,y,t) ).^.r* 0 rh 
*j=i 


From this we conclude that there exist B G L‘^{S, T*SQT*S) and (C, (p, g) G 
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L‘^{S, D{Uo)) such that 


symG'* ^ B + tSl^+Uo{C,ip,g), 


where 

B = B^ + {{MB%p + {MB^)p^)T'^QTf^ 

c = w + {{Mk)i,{Mk)2Y + {(j)i,(j)2y 

( \{Pi~kli + di) \ / |(Mi?^)3i \ ^ / dy^{Mk)2, \ 

5' = ( \{P2 + h + d2) ) I 1 2 ( j 


□ 

Remark 3.8. We were able to completely identify oscillatory part of two- 
scale limit of strain without using the convexity of the shell. The weak limit 
contained in the field B remains unknown. We will use the convexity in the 
construction of the recovery sequence by showing that in that case an arbitrary 
matrix field B can be recovered. In the plate case this is not completely true, 
but it is true that some appropriate class of the fields B can be recovered. 
However, it is commented that the class contains the minimizers that we 
need for the relaxation. 

We will just give the sketch of the proof of the lower bound since it follows 
standard arguments (truncation, Taylor expansion and lower semicontinuity 
of integral functional with respect to two-scale convergence). 

Proof of Theorem 13.41 - lower bound 

Without loss of generality we may assume that y^ dx = 0 and lim sup/^^g 
< oo. We only consider the case 7 G (0, cxd). The proof in other 
cases is similar. By using Theorem l3.5l we conclude that u G R^o(S') and that 
stronger convergence, given in that theorem are satisfied. Take an arbitrary 
u (£ u with boundary. By using Proposition 13.71 we have that there 
exists B G L‘^{S;T*S QT*S) and (Cw) ^ df{S,D{fA.f)) such that 

symG"^R + tS::+W,(C,p); 

As explained in |Neul2j (cf. [F,TMn2j for the corresponding argument in the 
homogeneous case), a careful Taylor expansion of lT(a; -|- tn{x),r{x)/e, I -\- 
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hG^{x)) combined with the lower semi-continuity of convex integral func¬ 
tionals with respect to weak two-scale convergence (see e.g. |Vis07( Proposi¬ 
tion 1.3]) and using relations flHOjl . flMjl yields the lower bound 


liminf = liminf 

h — ^0 h — ^0 



Q{x,t,y,B + tSl^{x) +U^{C,p)) dtdy dvols- 


si JIxY 


Minimization over B e L^(S', and (C,p) ^ L‘^{S, D{U^)) yields 


liminfh > [ Qy{x,Sli^{x))dvols■ 

h^O Jgi 

By using the properties of Qy (see Lemma 13.9p and exhausting u with the 
sequence {un)n such that for each n G N we have Wn d (u and Un has 
boundary we have the claim. 


3.4 Upper bound 

The proof of upper bound has the analogy with the planar cases. The crucial 
ingredient here is, as already commented in Remark 13.81 that every matrix 
held B can be recovered. It is easy to see, by using Korn’s inequality, that 
each hber of Ly{I x y) for 7 G [ 0 , 00] is a closed subspace of L^(/ x (V, 

Also by Korn’s inequality it is easy to see (see also [NeulOl [Neul2] i that the 
following coercivity bound is satished: 

||Wo(C,^)|li. > + \\g\\h) for all (C,^) G D{Uo), 

where the constant C depends on the embedding Analogous bounds are 
satished by Uy, for 7 G ( 0 , cxo] with the obvious norms on their respective 
domains of dehnition. 

The following lemma is analogous to [NV131 Lemma 2.10, 2.11]. 

Lemma 3.9. For 7 G [ 0 , cxd] we have the following 

(a) the function Qy : S x TS<^TS —>■ [0, cx)) is continuous, and there exists a 
constant C > 0 depending only on the energy density W and the surface 
S such that, for all x ^ S and all q E TfS Q TfS, 

C-^\q?<Q,{x,q)<C\q\\ (48) 

For fixed x E S, the function Qy{x) : TfS QTfS^ [0, cxd) is a quadratic 
form. 
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(b) there exists a bounded linear operator 

: L2(5, T*S ® T*S) L\S, T*S O T*S) x L\S, L^{I x 3^)), 
such that for almost every x E S we have 

Q^ix,q) = jj Q{x + tn{x),y,Il\[q]{x) + tq{x)+ T\f^[q]{x,t,y)^ dydt, 

whereYil^[q], n^[g] are projections of I\.^[q] on L‘^{S,TSqTS) and L‘^{S, Ly{Ix 
3^) respectively. Moreover, ifq G C^{S,TS®TS) thenll^[q] E C^{S,T*SQ 
T*S) X C^{S,L^iI xy)). 

Proof. We will just give the sketch of the proof for 7 G (0, 00 ). We introduce 
the space Vx = T*S Q T*S + x y) (the sum is orthogonal in the 
norm and thus the projections on the components are well defined and linear). 
Notice that for fixed x E S and q E TfS®TfS we can look the minimization 
equation (1551) as projecting the function tsymg on the linear space Vx in 
the norm determined by the quadratic form Q (which is equivalent to the 
norm). From this we obtain the uniqueness and the linearity of the 
minimizer for q G TfSQTfS. (a) can be now easily proved (see also |NV13] i. 

Also the part (b) is easily proved, hrstly for q G C^{S,TS ®TS) and then 
for q G L‘^{S,TS ® TS) by extension, using the inequalities in (1551) . The 
continuity of for q G C^{S,TS ® TS) follows from the property (Ql) 
stated in Lemma 13.21 □ 

The proof of the following lemma is given in [LMPini Lemma 5.6]. It uses the 
results in |Nir53j . We just augment it with the regularity statement that can 
be easily seen by the arguments given there and standard elliptic regularity. 

Lemma 3.10. Let a; C 6e simply connected with boundary. Let s : 
w —)■ 6e of class C^’^(cJ,M^). Assume additionally that s{oj) defines the 

local surface such the Weingarten map is strictly positive (or negative) 
definite up to the boundary, i.e., there exists C > 0 such that 

\/z' E u),\/v E < 85 ( 2 :')'^ ' ^ C!\v\‘^, 

(_x 

where (Ss )„;3 = dc^n^dys, for a, (3 = 1,2; Then for every 

B E L^(a;,Mgyj^) the system 

ql, = B, {ql,)ap = lidasdgw + dawdys), for a,(3 = 1, 2, 

has a solution w E If we assume that s G C^’"(a;,M^), B E 

C^’"(a;, Mgyj^), then we have that w E ^^’"(a;, M^). 
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Remark 3.11. The system ql^ = B arises naturally in the geometrical ques¬ 
tion of the problem of continuation of the infinitesimal bendings (see, e.g., 
IHLPl,^ \LMPlTf ). In \Nir5tff this system arises as one step linearization of 
solving the Weyl’s problem (finding a convex surface with prescribed metric) 
by a cotinuation method. 

Proof of Theorem 13.41 - upper bound 

• case 7 G (0, cx)): Without loss of generality we may assume that u is 
star-shaped with respect to some ball B{0,r) centered at the origin. 
Notice that, as a consequence of this, we have that for := Xu, (0 < 

A < 1), it holds u^ d u. Take an arbitrary u G By Remark 

I2.12l we have that u G M^). Dehne u G hju = uor. 

Let R G (7^’"(S', SO(3)) be such that RTs = Vu. Take arbitrary B G 
CII(S,T*SqT*S), (C,ple C((S,CH^ X X C((S,CH^ X 3;)). 

Dehne u G C'3’“(S,M3), b G R G S0(3)), (C,p) G 

C^(S,C\I X 3^,R2)) X C^(S,C\I X y)) by 

u = u o r, Bij = {B o ffjTi ■ Tj, R = Ro r. 

For hxed 0 < Aq < 1, sufficiently close to 1, dehne G C'^’"(AoCi;, R^) 
as the solution of the system = B\xou). For 0 < A < Aq, dehne 
u^,w^ G C'3'“(cJ,R3), B^ g 00(0;, R2yX^), R^ G ^^’“(o;, SO(3)) by 

u^(z') = j^u{Xz'),w^{z') = j^w{Xz'),B^{z') = B{Xz'),R^{z') = R{Xz'). 

In an obvious way, by pushing forward, dehne u^,w^ G 0^’"(S,R^), 
B^ G CK{S,T*SqT*S), R^ G C'1’“(R,S0(3)). Notice that 

\\u^ - = 0, hm \\B^ - R||i2(^_R2x^2) = 0, 

A,Ao—5-1 A,Aq—>^1 

lim ||R^ — R||_H-i(a;,so(3)) = 0, (49) 

A<Ao 
A, Aq—^1 

i.e., 

liin ||m^ — 'u||// 2(5 ^r 3) = 0, liin \\B^ — B\\l2(^s,t*sqt*s) = 0, 

A<Aq A<Aq 

A,Aq—^1 A,Aq—^1 

lim \\R^ — R||//1(5 ,so( 3)) = 0. (50) 

A<Aq 
A , Aq —^ 1 
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We define the recovery sequence 


il^{x + tn{x)) = u^{x) + thR^{x)n{x) + hw^{x) 

—th‘^ {dr^w^{x) ■ R^{x)n{x)) r"(x) 

+heR^{x) {Ca{xR,r{x)/e)T°‘{x) + p{x,t,r{x)/e)n{x)). 


Notice that 

lim \\y^ - M^||/fi(s,R3) = 0. (51) 

h->-0 

After a short computation we obtain 

{R^YVny'^ = I + thSly.{x) + h{R^{x)yVw^{x)Ts (52) 

—h {dr^w^{x) ■ R^{x)n{x)) t°‘{x) ® n{x) 
+hU^{(, p){x,t,r{x)/e) + o{h), 

where lim^^o II^^IU°° = 0. Here we used ffT5|l . flT^ . fl26|l . fl36|l . 

Notice that 

sym (^{R^Y'Vw^Ts — {dr^w^{x) ■ R^{x)n{x)) t°'{x) ®n{x)) = B^{x). 
By using the objectivity property (W3), fldOjl we easily obtain from fl52|l 


lim 

h^O 


1 

7^ 


W {x + tn{x),r{x)/e, Vhy{x + tn{x))) 


-Q {x + tn{x),r{x)/e, S0A(a;) + B^{x) p){x,t,r{x)/e)) 


0 . 


From this we easily obtain 


lim = 

h —^0 



5 Jixy 


Q{x + tn{x),y, SOa(x) + B^{x) +U^{C, p){x,t,y)) dydtdvo\s{x). 


By using diagonalization argument and letting A, Aq —)■ 1, we can con¬ 
struct a sequence (y^) such that lim/i_j,o |||/^ — m||_h'1 (5 ,r3 ) = 0 and 


lim h-^riy^) = 
h^o ^ 



Q{x + tn{x),y, S'^{x) + B{x) + p){x,t,y)) dydtdvolsix). 


s Jixy 
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Next we take Bi G L‘^{S, T*SqT*S), (Cz, Pi) G L‘^{S, ViU^)) by defining 

Bi = u][s:,], u,{Ci,pi) = u^^[s:,]. 

By using approximation with the sequence (5") C Cq{S,T*S O T*S), 
(C”,p’^)n C C^{S,C\I X y,R^)) X C^{S,C\I X 3^) such that 

lim ll-B” —-Bz||l 2 (s,t*s'©t* 5) = 0) p”) “ (Co Pz)I|l 2 (s' d(w^)) = 0. 

n^oo n—>-oo v w/ 

and the diagonalization argument we easily obtain the recovery se¬ 
quence. 


• case j = oo the only difference is in the relaxation part of the recovery 
sequence. We will just give the recovery sequence. For ((,p,c) G 
Ci(S.Ci(/,C'i(VRt)) X Ci(S,Ci(I,C'(y))) X Ci(S,Ci(I,R^)) we 
define 


y^(x + tn{x)) = u^{x) + thR^{x)n{x) -|- hw^{x) 

—th"^ [dr^w^{x) ■ R^{x)n{x)) r"(a;) 

+heR^{x) (Ca(x, t, r{x)/e)T°'{x) + p{x, t, r{x)/e)) 


+h‘^R^{x) ' Ca{x,s) dsr"-|- / 03(0;, s) dsn 



y 


The results from the previous case can be easily repeated. 


case 7 = 0, ©2 < /i < For {(, p, g) G Cl{S, C\y, R^))xC'2(^, C^{y))x 
Cl{S,Cl{I X 3^,R3)) we define 

y^{x -I- tn{x)) = u^{x) + thR^{x)n{x) + hw^{x) 

+e‘^R^{x)n{x)p{x,r{x)/e) —theR^{x) {dy^p{x,r{x)/e)T° 
—th? {dr^w^{x) ■ R^{x)n{x)) t°‘{x) 

+heR^{x) ((a(x,r(x)/e)T°'(x)) 

+h^R^(x) [2 f ga{x, s,r{x)/e) ds t°‘{x) 


gsi^x, s,r{x)/e) ds n{x) 


In this case we have 

{R^yVhy^ = I+ thSl,{x) + h{R\x)YVw\x)Ts (53) 

—h [dr^w^{x) ■ R^{x)n{x)) t°‘{x) 0 n{x) 

+edap{x, r[x)/e)n{x) 0 r"(a;) — edap{x, r{x)/e)T°‘{x) 0 n{x) 
+hUo{C, <P, 9 ){x, t, r{x)/e) + o{h). 
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Although in the expansion we have the dangeorus term of order e ^ h, 
the symmetric part of this term vanishes (see |Vell5] for the planar 
case). This means that 


{RjiR^yWhy’^ = I + thSl,{x) + h{R\x)YVw\x)Ts (54) 

—h (dr^w^{x) ■ R^{x)n{xy) t°‘{x) 0 n{x) 
RhUoiC, V?, g)ix, f, r{x)/e) + o{h), 

where R^ = exp(A^) and 

= eda^fix, r{x)/e)n{x) 0 t°‘{x) — eda^fix, r{x) / e)T°‘{x) 0 n{x). 

Again we can repeat the results from the first case (by using the ob¬ 
jectivity property (W3) and fl30|) b 
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